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INHOMOGENEOUS REFINEMENT EQUATIONS WITH RANDOM 

AFFINE MAPS 

RAFAL KAPICA, JANUSZ MORAWIEC 


Abstract. Given a probability space (n,A, P), random variables L,M: V, — >■ K. and 
g G L^(]R) we obtain two characterizations of these / S L^(]R) which are solutions of 
the inhomogeneous refinement equation with a random affine map of the form f{x) = 
In l■^(w)l/(^(w)a; - M{uj))P{duj) + g{x). 


1. Introduction 

Discrete inhomogeneous refinement equations of the form 

fix) = ^ Cnf{kx -n)+ g{x) 

n&, 

have been used in [8] for a construction of multiwavelets from a ’’fractal equation”, in 
[27] for a construction of boundary wavelets and in |28] for the study of convergence of 
cascade algorithms. Discrete inhomogeneous refinement equations and their continuous 
counterparts of the form 

fix)= / fikx-y)dg{y) +g{x) 

Jr 

were studied in [a to provide a uniform and a complete characterization of the existence 
of their distributional solutions to cover all cases of interest. Some other motivations for 
the study of discrete as well as continuous inhomogeneous refinement equations can be 
found in [6lClC2l[l3lElEni[2ll[22lE3l[2ll[25lE9lEni[3l]. 

Motivation to write this paper comes from applications, because from that point of view 
all results on the existence of ’’good” solutions of homogeneous as well as on inhomogeneous 
refinement equations are very important (see |3| where it is showed how nonexistence of 
’’good” solutions of a refinement equation can lead to anomalous behavior of numerical 
methods for a construction of wavelets). Let us note that refinement equations always 
have plenty of ’’bad” solutions, even extremely strange (see e.g. mini ESI). 

A common extension, both the above mentioned inhomogeneous refinement equations, 
is the following inhomogeneous refinement equation 

(1) fix)= [ \L{uj)\f{L{uj)x-M{u))P{duj) +g{x) 

Jn 
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with a random affine map <^ 9 : hi x M —)■ R defined by 
(2) ip{x,u;) = L{oj)x — M{u;), 

where L,M: Q ^ M. are given random variables and {Q,A,P) is a complete probability 
space. For more details on eqnation ([T]) with o' = 0 we refer the reader to the snrvey paper 
Cl and the references therein. 

First of all observe that if 5 ': R —)■ R is an integrable fnnction and f,h: R —)■ R are two 
representatives of a fnnction from L^(R) snch that ([I]) holds for almost all a: G R, then h 
also satisfies ([T]) for almost all a: G R; details are inclnded in [18]. This observation allows 
ns to accept the following definition: We say that / G L^(R) is an integrable solntion 
of eqnation ([!]), if ([ID holds for almost all x G R (with respect to the one-dimensional 
Lebesgue measnre). 

Given g G L^(R) let Vg denote the set of all integrable solntions of ([ID. It is obvions that 
Id; is a linear snbspace of L^(R) if and only if 5 ^ = 0. Moreover, if 7 ^ 0, then Vg is nniqnely 
determined by the snbspace Vq and by an arbitrary particular integrable solution of ([ID, 
if exists any. It is also easy to see that equation ([ID has at most one integrable solutions 
(which is clearly nontrivial) if and only if Vq = { 0 }. 

In this paper we are interested in the set Vg . More precisely, we obtain two character¬ 
izations of these functions / G L^(R) which belong to the set Vg . Moreover, we also get 
results guaranteeing that the set Vg is nonempty or consists of at most one member. 


2. Characterization of Vg by Fourier transforms 

From now on we fix a complete probability space (G, A, P), random variables L, M: G —)■ 
R, and a nontrivial function g G L^(R). 

If P(L = 0) > 0, then equation ([ID has exactly one integrable solution (see [IHl Corollary 
3.4], cf. [ini Remark 5.1]), and hence the set Vg consists of exactly one member. In this 
paper it is assumed that 

F(L = 0) = 0. 

Here and throughout, the Fourier transform of a function / G L^(R) is defined by 
= /Rexp{Rx}/(t)df. 

We begin with rewriting equation ([ID in the language of the Fourier transform. 

Lemma 2.1. Assume f G L^(R). Then f & Vg if and only if 

for every x G R. 

Proof. If / G Id;, then using the Fourier transform to both sides of ([ID and applying the 
Fubini theorem we conclude that ([3D holds for every x G R. 
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On the other hand, if ([3]) holds for every x G M, then applying the Fnbini theorem we 
obtain 


exp{itx}f{t)dt = 



Miu) +1 
exp < IX ———— 
L{u) 


f{t)dtP{du) + g{x) 



exp{mx}|L(a;)|/(L(a;)M — M{u))duP{du) + 'g{x) 


exp{itx} 


- M{u))P{du) + g{t) 


dt 


for almost all x G M. Since i i—>- Jq \L{uj)\f{L{u)t — M{u))P{du) + g{t) is an LOfunction 
having the same Fourier transform as f, if follows that ([T]) holds for almost all x G M. □ 


Putting X = 0 in ([3]) we get 

?( 0 ) = 0 , 

which is a necessary condition for the existence of integrable solutions of (ITi) . 

Given a measurable function : O ^ M define a sequence (\I/n)neN of random variables, 
acting on with real values, putting 


= ^{Un) 

for all n G N and zu = G Next with every h G L^(M) we associate a 

sequence {Ih,n)nm of real functions dehned by 






Mkjzu) 

(zaj ^ ‘' Lj^ ) 


h 


X 


Lii^'cu') • • • 


P°^{dw) 


for every x G M. 

Lemma 2.2. Assume that f G Then f & Vg if and only if 


( 4 ) 


N-l 

/(x) = If,N{x) + IgA^) + 

n=l 


for all N eN and x G M. 

Proof. For = 1 formula dH) coincides with formula ([3]), if we adopt the convention that 
X]n=i = 0 for every x G M. Hence Lemma 12711 implies that f E Vg if and only if for 

every x G M formula (jl]) holds with = 1. Therefore, to finish the proof it is enough to 
show that formula (jl]) holds for all A^ > 2 and x G M assuming that f E Vg . 
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Fix G M, / G and assume that (jl]) holds for every a; G M. Then by the induction 
hypothesis and Lemma [ 2.11 we obtain 


N 


fix) = 


MAw) 


exp <; j ^ 


f 


X 


P°^{dw) 


N-l 


^Ig,nix)+gix) 


n=\ 


/ exp zx V — 


Mkiw) 


(uj) ■ ■■Lk{w) 


exp < IX 


M{uj) 


f 

N-l 


X 


P{du) +g 


Li{w) ■ ■ ■ L]^{w)L{oj 

X 


Li (zu) ■ ■ ■ L]\f(^zu) 


P^{dw) 


Li{w) ■ ■ ■ Ln{w)L{u) 

+ ^h,nix) +gix) 

n=l 

N 

= If,N+lix) + ^Ig,nix) +gix) 

71=1 

for every x G M. The proof is complete. 

We are now in a position to formulate our hrst result of this section. 

Theorem 2.3. Assume f G L^(R). If 

(5) P{M = 0) < 1 
and the series 

'Z Mnizu) 

( 6 ) Z^[Zj) / T / \ T ( \ 

LAw ) ■ ■ ■ Lniw ) 

converges absolutely for almost all xd G then f ^ Vg if and only if 

« OO 

(7) /(x) =/(O) / exp{ixZoo{w)}P°°{dzu) + IgA^) + gjx) 

n=l 

for every x G M. 

Proof. According to Theorem 2.1 from [9], we conclude that 

(8) lim j j = 0 

77—^00 _L/-^ yZU) ' ' ' -L/n \'^) 

for all X G M and almost all vo G 12°°. 

Assume f G Vg . Then by the Lebesgue dominated convergence theorem we obtain 


□ 


lim If,Nix) = / exp{ixZooixxj)}fiO)P°°iduj) 

N^oo ' 


in° 
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for every x G M. This jointly with Lemma [2^ shows that the series ^g,n{x) converges 

for every x G M and, in consequence, ([7]) holds for every x G M. 

To prove the converse implication assume that ([7]) holds for every x G M. Then 




/ 


X 


L{uj) 


P{du) 


. M(u) 
exp <; IX—^ 
L{u) 

X 


n=l 


Liu] 


/(O) 


+ 9 


exp < IX 


^oo (^) 

Liu) 


P°°idw) 


X 


Liu) 


Pidu) 


I ) ■ ^oixo) , . 

j(U) / exp < zx -;-^—- p IX 


' 0 °° 


^oo i'^) 


f [■ Moiw 
y / exp<^zx——— 
I Loiw) 


Loiw) Loizu) 
M^iw 


P°^idw) 

MJw) 
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X 


Loixu). ■ ■ Ljiizu) 


+ 9 


X 


Loiw) 


P°°idw) 


/(O) / exp{zxZoo(ro)}P°°(dru) + 


n=l 


fix) -gix) 


for every x G M. Finally, by Lemma [2.11 we obtain that f GVg. 


□ 


Condition ([5]) is satished in all important applications of inhomogeneous (and homoge¬ 
neous) rehnement equations. In the boundary case 

(9) P(M = 0) = 1. 
we have the following result. 

Theorem 2.4. Assume f G L^(M). If holds and 

( 10 ) / ^og\Liu)\Pidu) <+00, 

then f & Vg if and only if 

OO 

(11) /(a^) = /(O) + 

n=l 


for every x G M. 
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Proof. From ([9]) we conclude that series ([6|) is well defined and Z^oi'co) = 0 for almost all 
vj G By the Kolmogorov strong law of large numbers we have 


lim 

n^oo 



exp<^ lim - Y^log|Lfc(ro)| 

1 n^oo n 


k=\ 


exp < / log |L(a;)|P((ia;) 


for almost all w G Thus applying flTOl) we infer that ([ 8 ]) holds for all x G M and almost 
all XD G Q°°. 

Assume f ^ Vg . Then by the Lebesgue dominated convergence theorem we obtain 
lim^v^oo If,N{x) = /(O) for every x G M. This jointly with Lemma YI 72 \ shows that the series 
Ig,n{x) converges for every x G M and, in consequence, flTT]) holds for every x G M. 

The proof of the converse implication is similar to that of Theorem 12.31 with an evident 
modification. □ 


Note that under conditions (jHl) and JH the set Vg consist of at most one member; 
indeed, if fi, f 2 e Vg , then /i — /2 G L^(M) and by Theorem 12.41 we have 

(/i - f2){x) = /i(x) - /2(x) = /i( 0 ) - 72(0) 

for every x G M, hence /i(0) = /2(0), and in consequence /i = / 2 - 

If (1T0|) holds, then the series ([6]) converges absolutely for almost all zu G 12°° provided 
logmax{|L“^(a;)M(a;)|, < +cx 2 and 

(12) P{cL + M = c) < 1 for every c G M 

(see [in], cf. [1]). Observe that flT^ rules out (jH]); indeed putting c = 0 in flT^ we obtain 
(Ej). It is known that replacing condition flTU]) by 


(13) — 00 < / log |L(a;)|P((ici;) < 0, 

Jq 

series ([6|) diverges if flT^ holds (see |9]). However, in the next result, convergence of series 
(]6|) is not relevant. 

Theorem 2.5. Assume f G P^(M). If flT3|) holds, then f ^Vg if and only if 

00 

(14) f{x) = Ig^n{,x) + g{x) 

n=l 

for every x G M. 


Proof. Assume f ^ Vg . Condition (IT^ jointly with the Kolmogorov strong law of large 
numbers implies 


lim 

TV—>-oo 


X 


= +CX) 


Lii^zu) ■ ■ ■ 

for all X G M \ {0} and almost all zu G 0°°. This jointly with the Lebesgue dominated 
convergence theorem and the Lebesgue-Riemann lemma yields 





lim 

TV^OO 


/ 


X \ 

Li{zu) ■ ■■Ln{zu)J 




0 
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for every x G M \ {0}. Therefore, Lemma YI?2 \ shows that the series converges 

for every x G M \ {0} and, in consequence, ffTTD holds for every x G M \ {0}. By the 
continuity of / and 'g we infer that ffTT|) also holds for x = 0. 

The proof of the converse implication is similar to that of Theorem 12.31 with an evident 
modihcation. □ 

Note that under condition ffT3|) the set Vg consists of at most one member, and moreover, 
if / G Vg , then /(O) = 0. 

We end this section pointing out that in the critical case 

(15) [ log\L{u)\P{du) = 0 

series (j6]) diverges in the case where (IT^ holds (see [1], cf. i) or converges in the case 
where ([9]) holds. But, independent of the case, it is well known that under assumption (IT^ 
equation (JT]) has completely different behavior than under assumption (ITOD or ffT3|) . The 
next two examples show that in the critical case equation ([1]) can have lots of integrable 
solutions. 

Example 2.6. Put P{uj) — L{uj) — u) and M{ia) — 0 for every w e 12 — {—1,1}. Then 
(Hg and ([9|) hold, and moreover, equation o takes the form 

(16) fix) = ^f{-x)+ ^f{x) + g{x). 

It is easy to see, taking g to he even and h G L^(R) to be odd, that the function f = h + g 
satisfies m- 

Example 2.7. Put P{oj) = 1, T(a;) = uj and M(a;) = 1 — a; for every oj & Q = { — 1,1}. 
Then fll5p and m hold, and moreover, equation o takes the form 

(17) /(x) = ^/(-x + 2) + ^/(x)+^(x). 

It is easy to see, taking g to have the graph symmetric with respect to the point (1,0) and 
h G L^(M) to have the graph symmetric with respect to the vertical line x = 1, that the 
function f = h + g satisfies (HHi, 

The above examples show why we are not interested in the critical case in this paper. 
More details on equation ([1]) in the critical case can be found in [T9] . 

3. Characterization of Vg by Radon-Nikodym derivatives 

Throughout this section we assume that L is positive almost everywhere on fl; the case 
where L is negative almost everywhere on hi is also acceptable for our purpose, but it will 
not be considered in this paper. 

It is known that every function from the set Vq is determined, up to a multiplicative 
constant, by the Radon-Nikodym derivative of their integrals over (— cxd,x], x G M with 
respect to the one dimensional Lebesgue measure (see m Section 3.4]). It turns out that 
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such a property also have all elements from the set Vg . To be more precise let us dehne a 
function G: M —)■ M by 

G{x) = / g(t)dt. 

J —OO 

It is clear that the function G is bounded. 

We have the following counterpart of Proposition 2 from |15] . 


Proposition 3.1. Assume f G Then f ^ Vg if and only if the function F: M —)■ M 

given by 

(18) F{x) = r mdt 

J —OO 

satisfies 


(19) 

for every a: G M. 


F{x)= [ F{L{uj)x-M{u))P{duj) + G{x) 

Jn 


Proof. Assume f & Vg . By the Fubini theorem we obtain 


F{x) = / f{t)dt = 



L{u) f {L{u)t — M{u))dtP{du) + / g{t)dt 


n J —OO 



L(u})x—M{lo) 


f{y)dyP{du) + G{x) = / F{L{uj)x — M{uj))P{du) + G{x) 


J —OO 


for every x G M. 

Conversely, if a function F: R —)■ M given by fflSj) satishes flT^ for every x G R, then 
Fubini theorem yields 


/ f{t)dt = F(x) = 

J —OO 

for every x G R. Hence f & Vg . 


L{u) f {L{u)t — M{u))P{du) + g{t) 


dt 


□ 


Now we need to write a precise formula for iterates (in the sense of Kuczma and Baron 
0 ) of the random map (p: R x H —)■ R given by ([2]); the iterates are dehned as follows: 

(p^{x,w) = (p{x,u!i) and ip"'~^^{x,'iu) = (p((p”(x, ro), a;„+i) 

for all X G R and zu = {oji,oj 2 , ...) G 
An easy calculation shows that 

n n n 

= x'^Lkiw) Tj(ro) 
k=l i=l J=i+1 

for all n G N, X G R and w G 12°°. Then, dehning a random variable Zn : 12°° —)■ R by 

n n—i 

Zn{uj) = - ^ Mj(tu) Lj{w) 
i=l j=l 


(20) 
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and a transformation cr„: —)■ by cr„ (wi, a;2, ■ ■ ■) = ■ ■ ■, <^1, ■ ■ ■) 'we obtain 

n 

an{w)) = X Lfc(ro) + Znivj) 
k=l 

for all n G M, X G M and w G It can be inferred from [10] (see also M) that assnming 
firsj) and 

(21) / logmax{|M(ci;) I, l}P((ia;) < cx) 

Jn 

the sequence of iterates ((/?"'(x, cr„(tJ 7 )))„gN converges for all x G M and almost all zu G 
and moreover, the limit is independent of x G M. Since for all x G M and n G N random 
variables 9?”(x, •) and 9?”(x, (j„(-)) have the same distributions, it follows that for every 
X G M the sequence •))neN converges in distribution to a random variable 0 , which 

is independent of x G M. In particular, the sequence {Zn)nm dehned by fl20|) converges in 
distribution to a random variable (j). Denoting by $ the probability distribution function 
of (f> we can formulate the following result. 


Theorem 3.2. Assume that M is integrable and let 


(22) 

[ L{uj)P{duj) < 1. 


Jn 

If 


(23) 

[ 9 {t)Ht)dt = g{0) 
Jr 


and G is Lipschitzian, then equation (IT^ has a Lipschitzian solution F: M —)• R. 

Moreover, if there exists f G L^(R) such that (ITS]) holds for every x G R, then f E Vg . 

The proof of Theorem 13.21 is based on assertion (iii) of Corollary 4.1 from |3]; this is the 
real reason for our assumptions about the random variable L. For the convenience of the 
reader, we repeat this assertion as a lemma formulated for requirements of our needs. 

Lemma 3.3. Assume that M is integrable and let (|2^ holds. If G is Lipschitzian and 
there exists Xq G X such that 

lim [ G{ip"^{xo,vj))P°^{dzu) = 0, 

n^oo 

then equation flT^ has a Lipschitzian solution F: R —)■ R. 

Proof of Theorem \3.A . We hrst note that condition (|22|) implies condition (IT^ and inte- 
grability of M forces condition 

For every n G N denote by the probability distribution function of the random 
variable Zn dehned by fl20l) . Then for every f G R, being a point of continuity of <h, we 
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have lim„^oo= ^{t)- Hence lim„^oo for almost all t G M. In 

conseqnence, by the Lebesgne dominated convergence theorem and fl2^ we obtain 

lim [ g{t)^n{'t)dt = [ g{t)^{t)dt = g{0). 

Jr Jr 

This jointly with the Fnbini theorem gives 


lim 

n—>-oo 


G{^^{0,w))P°°{dw) = lim / G{ip^{0,an{vj)))P°°{dvj) 


Q°° 




lim 

n^oo 


lim 

n—)-oo 


G{Zn{w))P°°{dw) 


g(t)dtP^ [dw) 


'n°° J-oo 


= lim 


g{t)dtP°°{dzu) 


g{t)P°°{d'aj)dt 



lim 

n^oo /TO /{roGO°°:t<Z„(ro)} 


lim / g{t)P°°{{zu G Q°° : t < Zn{w)})dt 

n^oo 

lim [ g{t)[l — P°°{{w G : Zn{w) < t})dt 

n^co 

g{0) - lim [ g{t)^n{t)dt = 0. 


Now the main part of the statement follows from Lemma 13.31 
The moreover statement follows from Proposition 13.11 


□ 


It is known that each absolntely continuous function F: M —)■ M of bounded variation 
with hma;^_oo F{x) = 0 has a suitably understood Radon-Nikodym derivative (due to the 
Hahn decomposition), i.e., F is represented by ffT5]l with / G L^(M). However, the following 
example shows that not all absolutely continuous real functions are represented in such a 
way. 


Example 3.4. There exists a bounded and Lipschitzian function F: M —)• M for which does 
not exist a function f G L^(R) such that flTSD holds for every x G M. 

Proof. Dehne a function F: M —)■ R by the Riemann improper integrals 

/ X 

h{t)dt, 

-OO 

where h: R —)■ R is given by 

f ^ for X G U“o(-2^- 

for X G U“=o(-2ri-2,-2n-1], 

[ 0 for X G (0, +cxo). 
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It is easy to show that F is Lipschitzian and bounded. Suppose, to derive a contradiction, 
that there exists / G L^(M) such that flTS]) holds for every x E X. Since f{t)dt = 
h(t)dt for all x,y eM., we have h = f almost everywhere on M. Hence h E L^(M), a 
contradiction. □ 
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